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RELATION OF LAND USE/LAND COVER TO RESOURCE DEMANDS 


SPECIFIC EMPHASIS: RESIDENTIAL ENERGY DEMAND 


SPECIFIC TASK: INVESTIGATE PREDICTIVE (FORCASTING) MODELS 


SUB-TASKS: 1. ECONOMETRIC MODELS OF RESIDENTIAL ENERGY 

DEMAND 

2. EVALUATION OF DETERMINANTS OF RESIDENTIAL 
ENERGY DEMAND VIS A VIS DERIVABILITY FROM 
REMOTELY SENSED DATA 

3. DATA BASES - ORGANIZATION AND INTEGRATION 

«. RESIDENTIAL LAND USE AND REMOTE SENSING, i.e., 
LAND USE/LAND COVER CLASSIFICATION 

5. LAND USE/LAND COVER CHANGE DETECTION 

6. LAND USE/LAND COVER PREDICTIVE MODELLING 



ECONOMETRIC MODELS OF RESIDENTIAL ENERGY DEMAND 


PURPOSE: 

UTILITY: 

FORM: 


SIGNIFICANT 


ISOLATE THE MAJOR DETERMINANTS OF DEMAND FOR ENERGY 


PREDICT CONSEQUENCES OF (i) NATURAL CHANGES SUCH AS 
CLIMATIC FLUCTUATIONS, AND (11) DELIBERATE POLICY 
ACTIONS BY DECISION-MAKERS SUCH AS PRICING POLICIES 


REGRESSION FRAMEWORK - 

Y = a + b^X-| + bgXg, ...» +b n X n> 

WHERE Y (DEPENDENT VARIABLE) = ENERGY DEMAND, X. 
(INDEPENDENT VARIABLES) = FACTOR DETERMINING ENERGY 
DEMAND, n = NUMBER OF INDEPENDENT VARIABLES, AND 
a and b. ARE MODEL PARAMETERS. 


NOTE: (i) INDEPENDENT VARIABLES (X.) ARE KNOWN TO 

INFLUENCE SIGNIFICANTLY THE VALUE OF Y 

(ii) PARAMETER 'a 1 IS A CONSTANT OF PROPORTION- 
ALITY AND PARAMETER ' b ' INDICATES THE CHANGE 
IN Y COMMENSURATE WITH A UNIT CHANGE IN X, i.e., 
MEASURE OF ELASTICITY 

EXAMPLE: Y = 3.0 -0.7X-, + 0.3X 2> 

WHERE Y = ENERGY DEMAND, X, = PRICE OF 
ENERGY, AND X £ = HOUSEHOLD 1 INCOME 


INDEPENDENT VARIABLES: 

PRICE OF FUEL (ELASTICITY OF DEMAND) 

HOUSEHOLD INCOME (ELASTICITY OF DEMAND WITH REGARD TO INCOME) 

PRICE OF SUBSTITUTE FUELS (CROSS ELASTICITY WITH REGARD TO 

ALTERNATIVE FUELS) 

HOUSEHOLD SIZE 
CLIMATIC CHARACTERISTICS 


PRICE OF HOUSEHOLD APPLIANCES SUCH AS HEATING AND COOLING 
APPARATUS 



ENERGY DEMAND DETERMINANTS AND REMOTE SENSING 


DETERMINANTS: PHYSICAL VS. NON-PHYSICAL ENVIRONMENTAL ATTRIBUTES 

PHYSICAL ATTRIBUTES: THEY ARE SCALE-DEPENDENT. ANDERSON LEVELS I AND II 

CAN BE DERIVED FROM LANDSAT DATA. 

URBAN LAND USE/LAND RESIDENTIAL - HIGH DENSITY 
COVER ATTRIBUTES: 

- MEDIUM/LOW DENSITY 

OPEN SPACE 

COMMERCIAL/ INDISTRIAL/ INSTITUTIONAL 
TRANSPORTATION - HIGHWAYS 

- RAILROADS 

- AIRPORTS 


THIS IS INADEQUATE. IT IS NECESSARY TO HAVE ADDITIONAL COLLATERAL INFORMATION. 
SUCH INFORMATION (SOCIO-ECONOMIC AND DEMOGRAPHIC ATTRIBUTES OF A POPULATION) 

IS OFTEN IN TABULAR FORM. IMAGE, MAP, AND TABULAR DATASETS MUST BE INTEGRATED 
INTO A GEO-BASED INFORMATION SYSTEM. 



Digital 

Codes 

First-Order Land Use/Land Cover 
Second -Order 

Third-Order 

X 

Urban and built-up land 

11 

Residential 

12 

Commercial and services 

121 

Recreational 

13 

Industrial 

14 

Extractive 

15 

Transportation, communications, and utilities 

151 

Utilities 

16 

Institutional 

17 

Strip and clustered development 

18 

Mixed urban 

19 

Open and other urban 

191 

Solid-waste dump 

192 

Cemetery 

2 

Agricultural land 

21 

Cropland and pasture 

211 

Nonirrigared cropland 

212 

Irrigated cropland 

213 

Pasture 

22 

Orchards, groves., and other horticultural areas 

23 

Feeding operations 

24 

Other agricultural land 

3 

Rangeland 

31 

Crass 

32 

Savannas 

33 

Chaparral (taken as brushtand) 

34 

Desert shrub 

4 

Forest land 

41 

Deciduous 

411 ! 

Deciduous/intermittent crown 

42 

Evergreen (coniferous and other) 

421 

Coniferous/solid crown 

422 

Co niferous/intcr mitten t crown 

43 

Mixed forest land 

5 

Water 

51 

Streams and waterways 

52 

Lakes 

53 

Reservoirs 

54 

Bays and estuaries 

55 

Other water 

6 

Nonforested wetland 

61 

Vegetated 

62 

, Bare 

7 

Barren land 

71 

Salt flats 

72 

Beaches 

73 

Sand other than beaches 

74 

Bare exposed rock 

741 

Hillslopes 

75 

Other banen land 

8 

Tundra 

81 

Tundra 

9 

Permanent snow and icefields 

91 

Pcrmanenc snow and icefields 




MINNEAPOLIS 1 

HOUSTON 2 

MILWAUKEE 3 

LOS ANGELES 4 

ORLANDO 5 












Commercial 

Core 

Industrial 

Core 

Commercial/ 

Industrial 

Strip 

Commercial/ 

Industrial/ 

Transporta- 

tion 

Commerce/ 

Industry 

Commercial/ 

Industrial/ 

Institu- 

tional 

Commercial 

Industrial/ 
New Con- 
struction 

High Density 
Single Fa- 
mily Res. 

Low Density 
Single Fa- 
mily Res. 

Mixed Single 
Multiple 
Family Res. 

Residential 

Residential 

(New) 

Mixed 

Residential 

Inner City 

Wooded 

Suburbs 

New 

Suburbs 

Other 

Suburbs 

Med . /High 
Density 
Residential 
Low Density 
Residential 

Residential 

Wooded 

Residential 

Urban Open 
Extractive 

| 


Undeveloped 

Urban 

Green Space 
Flood 

Channels & 
Extractive 

Undeveloped 


Woody Veg. 
Non Woody 
Veg. 

Trees 

Grassy 

Rural 

Chapparal 

Grassland 

Agricultural 

Trees 

Marsh 


Water 

Water 

Water 

Water 


L 


^Source ; 
^Source : 
^Source : 
^Source 
^Source 


Brown and Sizer, 1973 
Dornbach and McKain, 1973 
Mausel, Todd, and Baumgardner, 1974 
Bryant, 1976 

Hannah, Thomas, and Esparza, 1975 






















URBAN ACTIVITIES IDENTIFIABLE AT FOUR LEVELS OF INTERPRETATION 


ERTS-1 Satellite 
Imagery 

High Altitude Photography High Altitude Photography 
RB-57, 1:120,000 RB-57, 1:60,000 

Medium Altitude Photography 
Black-and-white, 1:15,R40 

Core 

Individual Structures 

Single Family Residential 

Housing Types 

Residential/Commcrcial 

Residential Areas 

Swimming Pools 
Apartment Complex 
Mobile Home Park 

High Rise Structures 
Garden Apartments 


Shopping Plaza 

Mobile Home Sales 

Pleasure Boat Sales 


Commercial Cluster 
Strip Commercial 
Administrative Buildings 
Schools 

University Complex 
Cemetery 

Parking Lots with Cars 

Building Under Construction 
Institutional Buildings 


Golf Course 

Boat Dock with Small 

Power Boat — Wake 


Baseball Diamond 
Drive-in Theater 
Marinas 

Boats 

Park 


Heavy Industry 

Junk Yard 

Power Plant — 


Tank Farm 

Extracting Industry 

Coal Piles 

Excavations 

Airports 

Highways 

Light Industry 

Excavating Industry 
Airport 

Terminal Building 
Aircraft Hangars 

Highway Interchanges 
Divided Highways 
Bridges 
Rest Areas 

Fabricating 
Processing 
Gas Storage 

Overhead Crane 
Water Pipes 
Open Storage Area 

Railroads 

Utilities 

R.R. Switching Yards 

Power Line Right of Way 
Secondary Roads 
Tertiary Roads 
Port Facilities — Ships 


R.R. Box Cars 


Source: Gary K. Higgs and M. Sullivan, ~A Compariiive Analysis of Remote Sensing 
Scak/System Attributes for a Multi-Level Land Use Classification System," 
Proceedings of the American Society Fall Convention (Falls Church, Va., 1973), pp. 
335-367. 
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DEFINITION: 


FORMAT: 


GEO-BASED INFORMATION SYSTEMS 


INTEGRATION OF DATASETS HAVING AREAL COLLECTION UNITS OF 
VARIED SIZE, SHAPE, AND FORMAT. 

LANDSAT IMAGE - REGULAR CELL-BASED (PIXEL) 

M.S.S. BANDS 4, 5, 6, AND 7 AND BAND RATIOS 

MAP FORMAT - POINT DIMENSION, e.g., ELEVATION, RAINFALL, INTERSECTION 
LINE DIMENSION, e.g., WATERCOURSE, RAOD, TRANSMISSION 

LINE 

AREA DIMENSION, e.g., CENSUS TRACT, NEIGHBORHOOD, 

LAND USE/LAND COVER CATEGORY 

TABULAR FORMAT - IRREGULAR AREAL UNIT FOR WHICH SOCIO-ECONOMIC 
INFORMATION IS AVAILABLE, e.g., AVERAGE FAMILY 
SIZE, AVERAGE FAMILY INCOME, AND AVERAGE NUMBER 
OF AUTOS PER FAMILY 


EXAMPLE: 


I.B.I.S. (IMAGE BASED INFORMATION SYSTEM) 




QUADRANGLE MAPS 




DATA 


IBIS 


DATA BASE 




BASED DATA 


Conceptualized Formation of an IBIS Data Base. 
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ISSUES 


LAND USE/LAND COVER CLASSIFICATION 


ACCURATE CLASSIFICATIONS OF REMOTELY SENSED DATA IS CENTRAL **** 
TO ALL APPLICATIONS TO WHICH THESE DATA ARE PUT. 

MANY ANALYSES AND MODELS ADD THEIR OWN VARIANCE TO THE DATA **** 
AND, THUS, IT IS CRITICAL TO BEGIN WITH HIGHLY ACCURATE DATA 
THAT CAN WITHSTAND SOME DEGRADATION DURING THE DERIVATION 
OF FINAL PRODUCTS. 

ASSOCIATED WITH THE ATTAINMENT OF ACCEPTABLE LEVELS OF ACCURACY: 

1. CLASSIFICATION SCHEME USED (DESIGNED FOR USE WITH REMOTELY SENSED 
DATA OR NOT). 

2. SPECTRAL AND SPATIAL FEATURE SELECTION AND POSSIBLE NEED TO 
REDUCE DIMENSIONALITY. 

3. IMAGE RADIOMETRIC (SUN ANGLE AND HAZE) AND GEOMETRIC (REGIS- 
TRATION) CHARACTERISTICS OF RAW AND PROCESSED REMOTELY SENSED 
DATA. GEOMETRIC TRANSFORMATIONS INVOLVE RE-SAMPLING. 

4. TEMPORAL SEQUENCE OF COVERAGE. 

5. GENERATION OF CLASS TRAINING STATISTICS AND AMPLING TECHNIQUES. 

6. OPTIMAL USE OF COLLATERAL ' DATA SUCH AS TERRAIN DATA WHICH CAN BE 


TREATED AS AN ADDITIONAL CHANNEL. 
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LAND USE/LAND COVER CHANGE DETECTION 


AIM: 


1. IDENTIFY THOSE PICTURE ELEMENTS THAT HAVE CHANGED OVER TIME 

2. A. DISTINGUISH CHANGE THAT IS OF INTEREST FROM THAT WHICH IS NOT 

B. CLASSIFICATION OF THE CHANGES OF INTEREST 


METHODOLOGY: 1. 

2 . 

3. 

4. 

5. 


IMAGE DIFFERENCING 
IMAGE RATIOING 
CLASSIFICATION COMPARISON 
PREPROCESSING PRIOR TO CHANGE DETECTION 
CHANGE VECTOR ANALYSIS 


1. IMAGE DIFFERENCING. 

A. PRECISELY REGISTER SCENES OF SAME AREA FOR TWO DATES 

B. SUBTRACT ONE IMAGE FROM THE OTHER ON A PIXEL-BY-PIXEL BASIS 

C. CREATE FREQUENCY DISTRIBUTIONS OF RADIANCE CHANGE FOR EACH BAND 

D. LOCATE THRESHOLD BOUNDARIES ON DISTRIBUTION TO SEPARATE CHANGE 
AND NO-CHANGE PIXELS (OFTEN DONE INTERACTIVELY' BY AN INTERPRETER 
FAMILIAR WITH THE AREA) 

2. IMAGE RATIOING. 

A. RATIO TRANSFORMATIONS TEND TO REMAIN INVARIANT UNDER VARYING 
CONDITIONS SUCH AS SHADOW, SUN ANGLE, AND SEASONAL REFLECTANCE 
DIFFERENCES 

B. NORMALIZED RATIO VALUES ON A PIXEL-BY-PIXEL BASIS IN ONE BAND 
ARE REPRESENTED BY A FREQUENCY DISTRIBUTION 

C. THRESHOLD SELECTION TAKES PLACE. RATIO VALUES SIGNIFICANTLY 
DIFFERENT FROM 1.0 ARE CONSIDERED CHANGE PIXELS 



IMAGE DIFFERENCING 
CHANGE DETECTION METHOD 
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3. CLASSIFICATION COMPARISON. 

A. SUCCESS DEPENDS SIGNIFICANTLY UPON ACCURATE CLASSIFICATION OF 
LAND USE/LAND COVER 

B. CLASSIFICATION OF MULTI-DATE, MULTI-SPECTRAL DATA WITH HOPE 
THAT A "CHANGE" CATEGORY WILL EMERGE 

C. LAYERED CLASSIFICATION REQUIRES MULTI-STAGE DECISION LOGIC. 
CLASSIFICATION STRATEGY IS BEST PERCEIVED AS A TREE DIAGRAM 

D. CLUSTERING COMPARISON. UNSUPERVISED CLASSIFICATION OF SCORES 
FOR TWO DATES AND COMPARION OF THE RESULTING GROUPINGS 

4. PREPROCESSING PRIOR TO CHANGE DETECTION. 

A. LOW FREQUENCY FILTERING - IMAGE SMOOTHING TO ENHANCE AREAS OF 
HOMOGENEITY AT EXPENSE OF HIGH FREQUENCY DETAIL 

B. HIGH FREQUENCY FILTERING - THIS ENHANCEMENT TECHNIQUE PRODUCES 
A SHARP VISUAL IMAGE BUT CONTAINS MORE "NOISE". GREATER EDGE 
ENHANCEMENT HIGHLIGHTING DISCONTINUITIES IN THE DATA 

C. PRINCIPAL COMPONENT ANALYSIS - REDUCES THE DIMENSIONALITY OF 
THE DATA BY COLLAPSING A SET OF CORRELATED VARIABLES INTO A 
SMALLER SET OF ORTHOGONAL VARIABLES 

5. CHANGE VECTOR ANALYSIS. 

A. WHEN LAND CHANGES ITS SPECTRAL APPEARANCE MAY CHANGE. PLOT 
THE VALUE OF TWO SPECTRAL VARIABLES (BANDS) FOR AN AREA BEFORE 
AND AFTER CHANGE ON A GRAPH 

B. VECTOR HAS BOTH MAGNITUDE AND DIRECTION: 

(i) MAGNITUDE INDICATES AMOUNT OF CHANGE; 

(ii) DIRECTION INDICATES TYPE OF CHANGE 

C. POSSIBILITY OF USING STANDARD DEVIATIONAL ELLIPSE ANALYSIS 




SPECTRAL VARIABLE X 



threshold 


threshold 


Decision 

T>*»ESHOLD 
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LAND USE/LAND COVER PREDICTIVE MODELLING 

LAND USE PATTERNS ARE A RESULT OF MISGUIDED SUBSIDIES, **** 
INSTITUTIONAL STRUCTURES, TECHNICAL CHANGE, CULTURAL LAG 
TIME AND, TO SOME EXTENT NON-QUANTIFIABLE INFLUENCES OF 
PUBLIC PREFERENCES. ALL THESE PARAMETERS INTERELATE TO 
PRODUCE A PATTERN. 

THE PATTERN IS COMPRISED OF A SET OF STATES (LAND USE 
CATEGORIES) AND THEIR INTERRELATIONSHIPS (SPATIAL PROX- 
IMITIES OR RULES GOVERNING THEIR INTERCHANGABILITY) . 

CONSIDER MODELS THAT ATTEMPT TO DESCRIBE AND EXPLAIN THE 
LAND USE CHANGE PROCESS. CONSIDER IN TURN: 

A. A TYPOLOGY OF LAND USE CHANGE MODELS 

B. OPERATIONAL MODELS INCORPORATING REMOTELY SENSED DATA 

C. SUGGESTED MODEL EXTENSIONS 
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A. A TYPOLOGY OF LAND USE CHANGE MODELS. 


I. INDEPENDENT MODEL. 


LAND USE IN CELL g.. AT TIME t+dt IS IN NO WAY RELATED TO THE LAND USE 

* J 

IN CELL g., AT TIME t. 

* J 


II. DEPENDENT MODEL (MARKOV MODEL). 

LAND USE IN CELL g.. AT TIME t+dt DEPENDS ON PREVIOUS LAND USE IN CELL 

* J 


g.. AT TIME t. 

* J 


t+dt 

ij 


F(g‘j) 


III. HISTORICAL MODEL (TIME SERIES OR LAGGED VARIABLE MODEL) 

LAND USE IN CELL g,, AT TIME t+dt DEPENDES ON THE SEVERAL PREVIOUS LAND 

* J 

USES IN CELL g... 

* J 

t+dt . t t-dt _t-2dt „t-kdt, 

9 ij • F( 9ij* 9lj • 3,'j ■ ■ • • • 9,j ) 

IV. MULTIVARIATE MODEL (MULTIPLE LINEAR REGRESSION AND DISCRIMINANT ANALYSIS) 
LAND USE IN CELL g., IS DEPENDENT ON SEVERAL OTHER VARIABLES ALSO AT -CELL 

* J 


ij 


t+dt rf t t t t t 

g.. = F(u-., v • • » w z-.) 

' j' ij’ ij 5 ’ 


V. GEOGRAPHICAL MODEL 

LAND USE IN CELL g., IS DEPENDENT ON LAND USE IN OTHER CELLS. 

* J 

t+dt _ p / t t 

9 ij 9 i±p» j+q 


A. EXTRAPOLATION-FILTERING MODEL 


- F(9 ?tpj+q> 

THIS MODEL CAN BE CHARACTERIZED BY A GEOGRAPHICAL QUIZ. 
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B. DYNAMIC GEOGRAPHICAL 


_t+dt _ p/_t n x 
q • « ~ h ( q . ..n. . ) 

lj*' 


WHERE n.. DENOTES ALL OF THE LAND USES IN THE NEIGHBORHOOD OF 

* J 


LOCATION i,j. 

THE MODEL HAS TV/O PARAMETERS: (i) - n - NEIGHBORHOOD, 

(ii) F - FUNCTION 


(i) NEIGHBORHOOD. THIS IS IMPORTANT BECAUSE IT DEFINES THE GEOG- 
RAPHICAL DOMAIN OF INFLUENCE. ONE COULD ASSUME THAT NEIGHBORHOOD 
DEFINITION VARIES BY CELL BUT IT IS EASIER AT ASSUME SPATIAL 
NEIGHBORHOOD STATIONARITY . ALTERNATIVELY, ONE COULD ALLOW 
NEIGHBORHOOD TO VARY IN SIZE, SHAPE, AND ORIENTATION AND BE A 
FUNCTION OF THE LOCATION OF THE CELL, i.e., p,q = F(i,j) 

(ii ) FUNCTION. ASSUME 5 LAND USE TYPES AND A CELL NEIGHBORHOOD OF 
5, i.e., (i,j), (i-1 ,j) , (i+l,j), (i,j-l), and (i,j+l). THUS, 
THERE ARE 5 STATES (S) AND 5 NEIGHBORS (N). 

A TRANSITION RULE SHOWS THAT ONE MUST CONSIDER S N = 3125 
CASES TO COVER ALL POSSIBILITIES. 

ASSUME SPATIAL ISOTROPY , i.e., POSITIONING OF NEIGHBORS DOES NOT 
COUNT AND SPATIAL STATIONARITY . THIS MEANS THAT THE SAME ENVIR- 
ONMENT (NEIGHBORHOOD) RESULTS IN THE SAME CONSEQUENCES OR THAT 
THE RULES DO NOT DEPEND ON WHERE YOU ARE. COMPARE CHESS RULES 
ARE PIECE-SPECIFIC BUT APPLY UNIFORMLY AT ANY LOCATION ON THE 


BOARD. 



MODEL I 


MODEL 13 


MODEL ID 


MODEL IY 


MODEL Y 






a 

□ 

□ 

□ 

B 

B 

□ 

□ 

a 

B 

B 

□! 

□ 

a 

B 

a 

a 

□ 

□ 

B 

B 

a 

B 

□ 

B 

B 

a 

a 

□ 

□ 

□ 

□ 

□ 

□ 

D 

a 


R R 

RAI *R Cl 

C C 


RJCRA *• C 


(2 R, II, \C,A) 
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B. OPERATIONAL MODELS INCORPORATING REMOTELY SENSED DATA 


I. LANDSCAPE MODELLING. 

DEFINITION: LANDSCAPE MODELLING ORGANIZES AND OVERLAYS INFORMATION FROM 
EXISTING MAPS, TABULAR SOURCES, AND THE ANALYSIS OF REMOTE 
SENSING IMAGERY INTO A COMPUTER FRAMEWORK. 


AIM: THE CURRENT THRUST OF LANDSCAPE MODELLING IS THE PROJECTION 

AND DISPLAY IN A MAP FORM OF THE FUTURE LANDSCAPE WHICH WOULD 
RESULT FROM THE CONTINUATION OF CURRENT LAND MANAGEMENT 
PRACTICES OR THE LACK THEREOF. 


METHODOLOGY:! . 

2 . 

3 . 

4 . 

5 . 

6 . 
7 . 


CELL-BY-CELL COMPARISON TO DETECT CHANGE 
MATRIX SHOWING FREQUENCY OF CHANGE BETWEEN CLASSES 
CONSTRUCTION OF TRANSITION PROBABILITY MATRIX 
MARKOV CHAIN MODELLING TO PREDICT FUTURE AGGREGATE CHANGE 
IDENTIFICATION OF DETERMINANTS OF CHANGE THROUGH DISCRIM- 
INANT ANALYSIS 

APPLICATION OF DISCRIMINANT MODEL TO ALL CELLS TO 

ESTIMATE NEXT MOST PROBABLE CHANGE 

ALLOCATE CHANGE BASED ON PRIORITY RANKING OF CELLS 



INPUTS: 



TOPOGRAPHIC MAPS 
VEGETATION MAPS 
SOILS MAPS 
GEOLOGIC MAPS 
TRANSPORTATION MAPS 


CENSUS TRACT MAPS OF 
POPULATION 
FAMILY COUNT 
HOUSING 
CAR OWNERSHIP 


SATELLITE IMAGES 
DIGITAL 
OPTICAL 

AIRCRAFT IMAGES 


PHOTOINTERPRETATION 
COMPUTER PATTERN CLASSIFICATION 

COMPUTES: 


SLOPE AND ASPECT PLANES 
SOCIOECONOMIC DENSITY PLANES 
MINIMUM-DISTANCE PLANES 
TRANSPORTATION ACCESS PLANES 
SOLAR-RADIATION PLANES 
LANDSAT TRANSFORMATIONS 


OUTPUTS: 

FUTURE LAND-USE SCENARIOS 
SITE DEVELOPMENT PLANS 
POWER PLANT SITE COMPARISONS 
ZONING ALTERNATIVES 

AND MANY OTHER APPLICATIONS TO 
ECONOMICS. HYDROLOGY. PLANNING, 
ENERGY. ETC. 


SIMPLE SCHEMATIC REPRESENTATION OF THE LANDSCAPE MODELING CONCEPT. 






7 LAND-USE SUBMODEL VARIABLES: 

1963 Photo Land -Use 

1970 Photo Land-Use 
1972 to 1973 US GS Photo Land-Use 
Two 1963 to 1970 Land-Use Changes 
Two 1963 to 1970 Alphanumeric Coded 
Land-Use Changes 

9 PHYSIOGRAPHIC SUBMODEL VARIABLES: 

— Topographic Elevation 
Topographic Slope 
Topographic Aspect 
USGS Surficia! Geology 

Solar Insolation for Land sat Imagery 
Four Landsat MSS -f Insolation Ratios 

5 TRANSPORTATION SUBMODEL VARIABLES: 

— Composite Minor-Road MD 
Composite Major-Road MD 
Freeway MD 

Freeway Interchange MD 
Built-Up Urban-Area MD 

17 SOCIOECONOMIC SUBMODEL VARIABLES: 
Four Population/Housing Densities 
Per Acre 

Five Population/Family/Housing- 
Unit Totals 

19£S Mean Family Income 
Median Housing-Unit Rent/Value 
One -/Two -/Three-Car Family Totals 
Total Census Tract Acreage 
Average Number of Cars Per Family 

10 LANDSAT IMAGE SUBMODEL VARIABLES: 

— MSS-4 (Visible Green) 

MSS-6 (Visible Red) 

MSS-6 (Solar Infrared) 

MSS -7 (Solar Infrared) 


CONCEPTUAL DIAGRAM OF THE LANDSCAPE MODEL OF THE DENVER STUDY AREA. 





Built-Up Urban Areas 


Minimum Distance to Built-Up Urban Areas 










CHANGES IN DENVER LAND USE, 1963 TO 1970. Entries show the number of 10 acre cells of each 1963 land use that 
converted to another land use by 1970. A blank denotes that no conversion occurred between the respective land uses during 
the period covered. Entries along the principal diagonal represent the number of cells of the land use that did not change to 
another land use between 1963 and 1970. 


Code* 11 


1963 Land-Use Type 


Residential 

Commercial and services 
Industrial 

Extractive 14 

Transportation, communl* 
cation, and utilities 

Institutional 16 

Strip and clustered . 

development 

Mixed urban 18 

Open and other 
urban land 

Cropland nnd pasture 21 

Orchards and other 22 

horticultural areas 4 

I'ccdinp operations 23 

Other agricultural land 
Deciduous forest land 
Streams and waterways 
Lakes 
Reserve rs 

Other water 

Vepclated nonforested 
wetland 

Sand other than beaches 
Other barren land 


1970 Land-Use totals 


11 6393 

\ 


12 13 


11 2 


16 17 


2 


Converting to 1970 Land Use 


21 22 23 24 


2 1089 


239 

52 

46 

35 

51 

42 

39 

3175 46 


1 

6 


u 








\ 






394 

27 

125 

139 

72 

31 

262 

351 14558 

2 

4 

22 

17 

5 


7050 I1R2 1061 636 865 3159 1655 0 3562 14624 


IH 101 641 175 









DENVER LAND-USE PROBABILITY TRANSITION MATRIX, 1963 TO 1970. Entries In the matrix denote the fraction of 
each 1963 land use that converted to another land use in 1970 as computed from table 4. A blank in the matrix denotes that 
no conversion occurred between the respective land uses during the period covered. A 1.0 in the matrix denotes that the area 
of that land use either did not change or increased by conversion to it from another land use between 1963 and 1970. Entries 
less than 0.001 are indicated by a dash. 


1963 bind- Use Type 


Residential 

Commercial and services 

Industrial 

Extractive 

Transportation, communi- 
cations and utilities 

Institutional 

Strip and clustered 
development 

Mixed urban 

Open and other 
urban land 

Cropland and pasture 

Orchards and other 
horticultural areas 

feeding operations 
Other agricultural land 
Deciduous forest land 
Streams and waterways 
Lakes 
Reservoirs 

Other water 

Vegetated nonforcstcd 
wetland 

Sand other than benches 
Other barren Lind 


B 










Converting to 1970 Lind Use 









11 

12 

13 

14 

15 

16 

17 

18 

19 2f 22 23 24 

41 

51 

52 

53 

55 

61 

73 

75 


1! 0.996 0.002 - 

\ 

12 0.002 0.988 

\ 

13 0.003 0.001 0.994 

\ 

14 0.006 0.011 0.959 

15 ^ 1.0 

16 0.004 0,001 

17 


0.995 


19 0.064 0.014 0.013 0.009 0.0 1 4 0.011 0.010 

21 0.025 0.002 0.009 0.004 0.002 0.016 

22 

23 

24 0.030 0.242 0.394 

41 

51 0.010 


O.OI6 0.141 


O.R49 0.012 

\ 

0.022 0.909 


0.061 0.091 


0.002 

o.noi o.ooi 


0,990 

\ 

0.998 


\ 

O.R I 3 






DENVER AGGREGATE LAND-USE PROJECTIONS BY A MARKOV TREND MODEL, 1963 TO 2068. Specific land-use 
total areas, such as those illustrated in figure 42, are detailed here. Photointerpreted second-order land-use changes between 
1963 and 1970 were used to drive this generalized trend model for the 24- by 24-mile Denver Metropolitan Area. The basic 
assumption of the Markov model is that the rates of changes are constant over time for all land-use classes. This model can 
project only aggregate land-use areas and cannot spatially predict the actual sites of conversion. 



Measured 

Acreages 







Projected Acreages 






1 

Land-Use lype 


1470 

1977 

1 984 

1991 

1998 

2005 

■g 

2019 

2026 

2033 

2040 

2047 

2054 

2061 


Residential 

m,:io 

70.500 

. 76.320 

.41,766- 

86.680 

91.280 

95,530 

99.450 

103.060 

106.400 

109.470 

112.310 

114.920 

117.330 

119.550 


Commercial and services 

11.020 

11.820 

. J1.JM 

j U.MNjf 
(•13, 760‘ 

13,960 

14.590 

15.180 

15.730 

16,260 

16.750 

17.210 

17.640 

18.040 

18.410 

1R.7R0 

19.120 

Industrial 

8.870 

10.610 

>12,250 

15.170 

16,470 

! 7,690 

18.810 

19.860 

20.830 

21.730 

22.570 

23.350 

24.070 

24.740 

25.160 

l.xtractlve . 

Transportation, contmuni- 

4/00 

6,360 

’ 7*800 

■ r . 

j> 9 AW 

|io,»6 

10.090 

10.980 

11,730 

12.350 

12.860 

13.260 

13.580 

13.810 

13,980 

14.080 

M.I30 

14.120 

cations, and utilities 

7.200 

8.650 

*20 

11,890 

12,810 

13,660 

14.450 

15.190 

15.870 

16.500 

17.090 

17.630 

18.140 

18.610 

19 040 

Institutional 
Strip and clustered 

J 1.250 

31.540 

31, M0 

'.'•1 • 

32,136 

' ■ C ■ 1 

32,330 

32.500 

32.620 

32.710 

32.770 

32.810 

32.810 

32,790 

32,740 

32.680 

32.600 

32.500 

development 

1.1,500 

16.550 

19,310 

2M50 

24,190 

26.3.10 

28.300 

30,110 

31.770 

33.300 

34,710 

36.010 

37.200 

38.300 

39. .120 

40 260 

Mixed urban 

40 

0 

o 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

n 

Open and other urban 

.17,4)0 

35.620 

33.800 

31.980 

30.200 

28.470 

26.RI0 

25,220 

23,700 

22.270 

20.920 

19.650 

18.460 

17.340 

16.310 

15 140 

Cropland and pasture 
Orchards and other 

1 60,040 

146.240 

133/>00 

122,120 

: 

111.670 

102.150 

93.500 

85.610 

78,440 

71.900 

65,940 

60.510 

55.560 

51.050 

46940 

4 3.180 

horticultural areas 

60 

60 


. -r.60 

60 

60 

60 

60 

60 

60 

60 

6(1 

60 

60 

60 

60 

l eedinp operations 

20 

20 

20 

i,‘ 20 

■ ■ l._ 

20 

20 

20 

20 

20 

20 

20 

20 

20 

20 

20 

20 

Other aprlciiltural land 

330 

70 

30 

70 

20 

20 

20 

i 

10 

10 

10 

10 

10 

10 

in 

10 

10 

Deciduous forest land 

urn 

180 

180 

/ 186 

180 


180 

180 

180 

180 

180 

180 

IKO 

180 

180 

180 

Streams and waterways 

460 

1.010 

1.060 

MOO 

1,140 

0 

1.210 

1.250 

1.280 

1.300 

1.330 

1.360 

1 ,J80 

1.400 

1.4 20 

1.440 

Lakes 

5,910 

6,410 

6,880 

7,340 

7.800 

8.260 

MOO 

9.140 

9.580 

10.000 

10,4 20 

10.830 

11.240 

1 1.630 

12.020 

1 2.400 

Reservoirs 

1.580 

J.750 

1,910 

2,050 

2.180 

2,300 

2.400 

2.500 

2.59() 

2.680 

2.750 

2.820 

2.890 

2.950 

3.000 

3. ft Ml 

Other water 
Vcpctaled nonforeslcd 

50 

50 

so 

50 

50 

50 

50 

50 

50 

50 

50 

50 

50 

50 

50 

50 

wetland 

1.710 

1.710 

1,7 10 

1,710 

1.710 

1.710 

1.710 

1.710 

1,710 

1.710 

1.710 

1.710 

1.710 

1.710 

1.710 

1.710 

Sand othet than beaches 

, 640 

520 

420 

340 

280 

230 

ISO 

150 

120 

100 

RO 

70 

50 

40 

30 

30 

Other barren land 

18.870 

18.920 

18.460 

19,000 

19.040 

19.070 

19,090 

19.110 

19.1.10 

19.140 

19.150 

19,160 

19.170 

19.170 

1 9.J70 

19.160 

Grand Total 
r— i _ 

568.640 

368.640 

368.640 

368,640 

368.640 

368.640 

3 68. 640 

368.640 

368.640 

368.640 

368.640 

368.640 

368.640 

368.640 

368.640 

168.640 


I I Oh«r*«lk>n interval. □ Region of possihle reality. 


































TIME (year) 

PREDICTION OF FUTURE TRENDS IN THE AMOUNT OF 
OPEN SPACE AND COMPETING LAND USE IN THE DENVER METROPOLI- 
TAN AREA. A constant matrix of transfers, P = | from land use, i f to 
land use, j, was assumed. Acreages displayed relate to the original area of 24- by 
24-statute miles = 576 square statute miles. 


2068 





PREDICTED 
1977 LAND- 
USE MAP 



PROPOSED COMBINATION MARKOV AND LINEAR DISCRIMINANT MODELS FOR IM- 
PROVED SPATIAL-CHANGE PREDICTION. The Markov trend model provides the correct number of 
change cells by type. These can be selected from a sorted list of discriminant-computed posterior probabili- 
ties of change. The selected cells can be assembled into a map of future land use. Spatially registered Land- 
sat digital imagery can serve as future land-use inputs in lieu of the 1963 to 1970 aerial photography. 
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SPATIAL ANALYSIS AND REGIONAL LAND USE PATTERNS. 

DEFINITION: THE TERMS SPATIAL ANALYSIS AND SYNTHESIS ARE DEFINED AS 

THE SEPARATION OF A LAND USE PATTERN INTO INDIVIDUAL 
COMPONENTS (SPATIAL PROXIMITIES) AND THE RECONSTRUCTING 
OF THESE COMPONENTS FOR USE IN PROJECTING LAND USE PATTERNS. 

AIM: NEED TO PROJECT FUTURE LAND USE PATTERNS. 


METHODOLOGY: 1. 

2 . 

3. 


4. 


5. 


BASED ON QUANTIFICATION OF A LAND USE PATTERN 
OPERATION DEFINITION OF SPATIAL RELATIONSHIPS 
DERIVATION OF A WEIGHTED INDEX BASED ON FREQUENCIES 
OF OCCURRENCE FOR EACH DISTANCE FROM ALL OTHER LAND 
USES 

IDENTIFICATION OF "TRIGGERING" AND "CONSTRAINING" 
FACTORS 

ALLOCATION PROCEDURE OF AGGREGATE DEMAND FOR LAND TO 
SPATIAL DOMAIN 


OPERATION: 


1. SELECTION OF 15 LAND USE/LAND COVER CLASSES CODED 
AS PRESENCE/ABSENCE ON CELL-BY-CELL BASIS 

2. CREATION OF SPATIAL PROXIMITIES MATRIX FOR EACH LAND 
USE/LAND COVER CATEGORY WITH 19 DISTANCE INTERVALS. 
ACCOMPANIED BY A MAP SHOWING ISO-PROXIMITY LINES 

3. TABLE OF PROXIMITY OCCURRENCES OF ANY LAND USE/LAND COVER 
CATEGORY TO ALL OTHER CATEGORIES BY DISTANCE INTERVALS 

4. PROJECTION TECHNIQUE. MODEL MUST DISCRIMINATE UNDEV- 
ELOPED AREAS HAVING THE SAME SPATIAL CHARACTERISTICS 
AS AREAS CURRENTLY OCCUPIED BY A GIVEN LAND USE. THUS, 
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A CELL WITHOUT URBAN DEVELOPMENT FOUND TO POSSESS THE 
LARGEST NUMBER OF SPATIAL ATTRIBUTES THAT MOST CLOSELY 
APPROXIMATE THE ATTRIBUTES OF DEVELOPED CELLS IS 
ASSUMED TO DEVELOP IN THE SAME MANNER. 

5. CROSS REFERENCE FILES SHOWING (i) CELLS CONTAINING A 
PARTICULAR LAND USE CATEGORY; (ii) CELLS HAVING HIGH 
WEIGHTING FOR SOME CATEGORY; (iii) LAND USE CATEGORY 
DEMAND (EXPRESSED IN ACRES OR NUMBER OF CELLS); (iv) 
CELLS CONTAINING CONSTRAINING ATTRIBUTES. 














TABLE OF FREQUENCIES COUNTS OF MEDIUM DENSITY 
HOUSING VERSUS DISTANCES* FROM OTHER LAND USES 


DISTANCES 


LAND USES 

0 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

ii 

12 

13 

14 

15 

16 

17 

18 

Water 

34 

90 

82 

77 

106 

66 

91 

74 

53 

43 

19 

16 

10 

13 

8 

0 

4 

2 

1 

3 j 

Institutions 

176 

379 

125 

45 

40 

20 

11 

B 

B 

0 

0 

2 

2 

B 

B 

0 

0 

0 

i 

0 

Low Density Housing 

m 





B 

5 

B 

5 

3 

0 

0 

0 

0 

0 

0 

0 

0 

i 

0 

High Density Housing 




93 

83 

26 

23 

14 

13 

Bil 

8 

6 

9 

3 

6 

•8 

12 

D 

8 

Recreation 

77 

186 




46 

24 

B 

11 

13 

B 

B 

3 

2 

6 

6 

B 

n 

6 

Railroads 

83 

184 

115 



63 

62 

51 

26 

m 

8 

6 

3 

B 

0 

0 

0 

0 

0 

Industry 

23 



97 1 

mm 

75 

76 

52 

46 

40 

17 

11 

3 

2 

B 

■ 2 

B 

]. 

0 

Commerce 

101 



89 


32 

20 

* 

> 

6 

B 

3 

B 

B 

B 

0 

0 



Rural (major) 

153 



B 

41 

23 

14. 

" 

w 

7 U/' 

14 


2 

0 

1 

0 

B 

0 

0 

0 

0 

Rural (minor) 



m 

69 


50 

B 


.35 

43 

1 

1 

30 

37 

31 

33 

24 

16 

9 

6 

0 

Controlled Access 

5 

B 

8 

5 

B 

15 

/ j 

14 

jJM 

B 

13 

8 

15 

1 

14 

19 

15 

B 


12 

’ 16 

Uncontrolled Access 

133 

240 

125 

95 

78 

29 

IB 

m 

17 

. 15 

9 

8 

6 

6 

B 

B 

2 

0 

1 

Transportation 

(interchanges) 

5 

26 


36 

H 

52 

50 

38 

47 

65 

51 

48 

36 

40 

27 

27 

21 

10 

ii 

Transportation (interstates) 

20 

55 

45 

36 

61 

53 

46 

39 

50 

63 

43 

38 

29 

29 

31 

20 

16 

12 

9 

Sewer Service Areas 

217 



72 

57 

28 

15 

15 

12 

14 

8 

11 

B 

5 

5 

B 

3 

5 

0 

Rivers 

130 




108 

31 

13 

B 


0 

0 

0 

B 


0 

0 

0 

0 

0 


•Distances equal 500 meters. 
CK/lmt 2/8/77 revised 11/23/77 
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C. SUGGESTED MODEL EXTENSION. 

MODEL SHOULD POSSESS TWO MAJOR ATTRIBUTES: 

I. ABILITY TO PREDICT THE MAGNITUDE OF AGGREGATE SYSTEM-WIDE 
LAND USE/LAND COVER CHANGE 

II. ABILITY TO PREDICT THE MOST LIKELY GEOGRAPHICAL LOCATION WHERE 
CHANGE WILL OCCUR. 


I. A. STOCHASTIC APPROACH. 

ATTENTION SHOULD BE PAID TO HOW WELL THE PHENOMENA BEING MODELLED 
APPROXIMATE THE MATHEMATICAL AND STATISTICAL REQUIREMENTS OF THE 
MODEL, e.g., STATIONARITY OF THE TRANSITION PROBABILITIES AND ORDER 
OF THE PROCESS. 

B. USE OF EXOGENOUS (TRIGGERING) VARIABLES WHICH ARE DERIVED FROM 
DEMOGRAPHIC/ECONOMIC MODELS DESIGNED TO PREDICT FUTURE LAND RE- 
QUIREMENTS. THIS APPROACH REQUIRES A DETERMINATION OF TRIGGERING 
FACTORS, CONSTRAINING FACTORS, AND SUITABILITY CRITERIA. 

II. A. MULTIVARIATE APPROACH WHICH DETERMINES THE RELATIONSHIPS BETWEEN 

CHANGE AND A NUMBER OF INDEPENDENT VARIABLES VIA MULTIPLE REGRESSION 
AND DISCRIMINANT ANALYSIS. 

B. USE OF LATERAL DEPENDENCIES AND THE CALCULATION OF CONDITIONAL 

PROBABILITIES ALONG THE LINES OF THE "GEOGRAPHICAL MODEL" DESCRIBED 
EARLIER. ATTENTION MUST BE PLACED ON THE DEFINITION OF THE CONCEPTS 
"NEIGHBORHOOD" AND "FUNCTION". 
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THE ABOVE CONCERN CAN BE FORMALLY EXPRESSED AS FOLLOWS-: 

P K = F ( pT K* P K)» 

WHERE pj IS THE TIME TRANSITON MATRIX DESCRIBING THE PROBABILITY 
THAT THE NEXT CELL STATE WILL BE K, p[! IS THE CONDITIONAL PROBABILITY 
THAT CELL TYPE K WILL OCCUR, GIVEN THE STATES OF THE NEIGHBORING 
CELLS, AND F IS A FUNCTION THAT COMBINES THE TWO PROBABILITIES. THE 
EXACT SPECIFICATION OF THE FUNCTION CAN VARY, e.g., 

F(pJ.pJ) = Pjj . OR 

F(P*.pg) = + W 2 P[f, WHERE 

W 1 AND W 2 ARE APPROPRIATE WEIGHTS. 






U.S. BUREAU OF THE CENSUS 
TABULAR INFORMATION FILES 


PUBLIC UTILITY 
COMPANY REVENUE 
AREAS 



PIXEL 


REMOTELY SENSED 
IMAGERY 


AUTOMATIC OR 
MANUAL CLASS- 
IFIER 





ATTRIBUTE LIST 
PER PIXEL: 


. CENSUS 01 
. CENSUS 02 


ENERGY 01 
ENERGY 02 
LANDUSE 01 
LANDUSE 02 
LANDUSE 03 




THEMATIC LAN DUSE MAP 



POPULATION 



POPULATION 





THEMATIC LANDUSE MAP 
WITH K CATEGORIES FOR 




LANDUSE CHANGE 
MODEL 

t 





ECONOMETRIC ENERGY 
DEMAND MODEL 


t 


SOCIOECONOMIC ATTRIBUTES 
FROM CENSUS OR ESTIMATES 


FUTURE ENERGY DEMAND: 

(A) AGGREGATE; 

(B) SPATIAL PATTERN 
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APPLICATIONS: 


MODEL ELEMENTS: 


MARKOV CHAIN MODELLING 


MODELLING MOVEMENT IN TERMS OF GEOGRAPHICAL MOVEMENT OR IN 
TERMS OF MOVEMENT FROM ONE "STATE" TO ANOTHER. A SYSTEM IS 
IN A PARTICULAR STATE AT ANY POINT IN TIME, i.e., HAS A 
PARTICULAR STRUCTURE OR SET OF RELATIONSHIPS AMONG THE PARTS. 

MARKOV CHAIN MODELS ARE NEAT AND ELEGANT CONCEPTUAL DEVICES 
FOR DESCRIBING AND ANALYZING THE NATURE OF CHANGE AND MAY BE 
USED TO FORCAST FUTURE CHANGE. 

PLANNERS AND ADMINISTRATORS CONCERNED WITH THE PRESENT AND 
FUTURE SPATIAL ORGANIZATION OF THE LANDSCAPE ARE INTERESTED 
IN THE PROBABILITIES OF CHANGE. THE MARKOV CHAIN MODEL IS 
A STOCHASTIC TYPE OF MODEL. 

1. TRANSITION PROBABILITY MATRIX WHERE ELEMENTS IN EACH ROW 
SUM TO 1.0. THE PRINCIPAL DIAGONAL OF THE SQUARE MATRIX 
INDICATES THE PROBABILITY OF NO CHANGE IN THE STATE. 

2. INITIAL STATE VECTOR DESCRIBES THE STRUCTURE OF THE ENTIRE 
SYSTEM AT ANY POINT IN TIME, e.g., PERCENT OF LAND IN 
EACH OF A NUMBER OF LAND USE/LAND COVER CLASSES. THE 
ELEMENTS OF THE VACTOR SUM TO 1.0. 

3. A CHANGE IN THE STATE VECTOR IS ACCOMPLISHED BY MULTI- 
PLYING THE INITIAL STATE VECTOR BY SUCCESSIVELY HIGHER 
POWERS OF THE TRANSITION PROBABILITY MATRIX. 

4. THE SIMPLE MODEL ASSUMES: (i) A CONSTANT POPULATION; 

(ii) A SET OF TRANSITION PROBABILITIES; (iii) THESE 
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PROBABILITIES REMAIN CONSTANT OR STATIONARY ; AND 
(iv) THE FIRST ORDER PROPERTY, i.e., THE SYSTEM 
CHANGING FROM A GIVEN STATE S 1 AT TIME t Q+1 DEPENDS 
ONLY ON THE STATE S i AT THE TIME t Q AND IS INDEPENDENT 
OF THE STATES OF THE SYSTEM PRIOR TO t Q . 

5. PROPERTIES OF REGULAR FINITE MARKOV CHAINS. 

(i) TRANSITION PROBABILITIES ARE NON-NEGATIVE AND THOSE 
IN EACH ROW SUM TO 1.0 AND THE TRANSITION MATRIX 

IS REGULAR IF FOR SOME POWER OF THE MATRIX THERE 
ARE ONLY POSITIVE NUMBERS. 

(ii ) CONCEPT OF EQUILIBRIUM IS EXPRESSED IN TWO THEOREMS: 

.a. IF P IS A TRANSITION MATRIX FOR A REGULAR MARKOV 

CHAIN THEN - 1 THE POWERS OF P APPROACH A MATRIX A 

2 EACH ROW OF A IS THE SAME PROBABIL- 
ITY VECTOR 'a' 

3 THE ELEMENTS OF 'a* ARE ALL POSITIVE 

.b. IF P IS A TRANSITION MATRIX FOR A REGULAR MARKOV 
CHAIN AND A AND 'a' ARE AS STATED ABOVE, THEN THE 
UNIQUE VECTOR 'a' IS THE UNIQUE PROBABILITY 
VECTOR SUCH THAT , a , .P= , a'. THE MATRIX A IS 
DEFINED AS THE LIMITING MATRIX . 

(iii) THE LIMITING MATRIX DESCRIBES THE AVERAGE STATE OF THE 
SYSTEM AND THE ASSOCIATED PROBABILITY VECTOR 'a' HOLDS 
THE SYSTEM IN EQUILIBRIUM. 

(iv) IN MARKOV CHAIN ANALYSIS, FOR MODELLING PURPOSES, THE 
EQUILIBRIUM DISTRIBUTION IS OF INTEREST NOT AS A 
FORCAST OF THE FUTURE STATE OF THE SYSTEM BUT AS A 



-18- 


PROJECTION OF WHAT WOULD BE IF THE OBSERVED PATTERN 
OF MOVEMENT CONTINUED UNHAMPERED. 

(v) THE FUNDAMENTAL MATRIX (Z) DESCRIBES HOW THE SYSTEM 
APPROACHES EQUILIBRIUM FROM AN INITIAL DISTRIBUTION. 
IT CAN BE USED TO CALCULATE A MATRIX OF MEAN FIRST 
PASSAGE TIMES , i.e., THE TIME IT TAKES ON AVERAGE 
TO MOVE FROM ONE STATE TO ANOTHER. 


QUESTIONS ABOUT THE APPLICATION OF THE MARKOV MODEL. 

1. IDENTIFICATION OF THE SPECIFIC ORDER PROPERTY SHOULD BE MADE. MOST 
APPLICATIONS ASSUME A FIRST ORDER. CRITERIA EXIST TO TEST SUCH AN 
ASSUMPTION. 

2. STATIONARITY OF THE PARAMETERS, i.e., THE ESTIMATED TRANSITION 
PROBABILITIES ARE FIXED OR CONSTANT THROUGHOUT THE PREDICTIVE PERIOD. 
STATISTICAL TESTS EXIST. 



Estimated probability transition matrix of commercial activities for downtown Denver. 
First row, 1947-1951; second row, 1951-1963; third row, 1963-1967; fourth row, 1967-1971. 


Transition from Transition to a 



1 

2 

3 

4 

5 

6 

7 

8 

9 

Financial, 

0-798 

_ 

_ 

0*125 

0-056 

_ 

0-014 

_ 

0 007 

real estate, 

0*724 

0 052 

0-013 

0-089 

0019 

0013 

0-038 

- 

0 051 

and insurance 

0-747 

- 

- 

0-070 

0-014 

0014 

0*042 

- 

0*112 


0 806 

- 

- 

0015 

0 029 

0-045 

0 059 

- 

0-045 

Professional 

0-026 

0-615 

_ 

0-218 

0-038 

_ 

0 027 

_ 

0-080 

services 

0-080 

0-400 

0-120 

0-040 

0-080 

0-040 

0160 

- 

o-oso 


- 

0-625 

- 

0-042 

0-208 

- 

0-125 

- 

- 


- 

1-000 

- 

- 

- 

- 

- 

- 

- 

Commercial 

- 

- 

0-919 

0-054 


_ 

0-027 

_ 

- 

residential 

0-095 

- 

0-667 

- 

- 

0 048 

0-190 

- 

- 


- 

- 

0-727 

- 

- 

- 

0*181 

- 

0 091 


- 

- 

0-682 

- 

- 

- 

0-046 

- 

-0-272 

General retail 

0-029 

0-003 

_ 

0-830 

0-024 

0-016 

0-012 

_ 

0 031 


0-033 

0-028 

0*002 

0-539 

0-040 

0 022 

0-191 

- 

0144 


0 045 

0-003 

• 

0-719 

0*010 

0-023 

0-069 

0-005 

0-126 


- 

0 003 

0 002 

0-703 

0-014 

0-028 

0-023 

0016 

0-217 

Personal 

0-030 

- 

- 

0146 

0-732 

0-021 

- 

_ 

0-098 

and commercial 

0-055 

- 

0-004 

0*164 

0-535 

0-008 

0*133 

- 

0102 

services 

0-057 

0-020 

- 

0-062 

0-647 

- 

0-095 

- 

0*119 


- 

- 

- 

0-072 

0-619 

0 041 

0 031 

- 

0-237 

Eating 

0-012 

- 

_ 

0*060 

0-036 

0-811 

0019 

. 

0 062 

and/or drinking 

0 069 

- 

0011 

0-042 

0-070 

0-643 

0168 

- 

0 049 


0-050 

- 

- 

0-033 

0-004 

0-767 

0 083 

- 

0 063 


0-009 

- 

- 

0-027 

0-009 

0-691 

0-036 

- 

0-227 

Parking 

0-067 

_ 

- 

0-022 

0011 

0-011 

0-845 

. 

0-044 


0-042 

- 

0-013 

0 042 

0-013 

0-030 

0-795 

0-013 

0-065 


0 033 

- 

- 

0-033 

- 

- 

0-900 

- 

0-033 


- 

- 

- 

- 

- 

- 

0-946 

- 

0-054 

Entertainment 

- 

_ 

- 

_ 


0-042 

_ 

0-958 

- 


0-037 

- 

- 

0-111 

- 

0-037 

0*185 

0-630 

- 


0-154 

- 

- 

- 

- 

- 

- 

0-846 

- 


- 

- 

• 

- 

- 

- 

0- 167 

0-666 

0167 

Other 

0-042 

0-004 

0-002 

0 081 

0-065 

0010 

0-007 

- 

0-798 


0-026 

0018 

0-018 

0-215 

0-044 

0 062 

0-106 

- 

0-499 


0-036 

0-001 

- 

0-066 

0-005 

0-009 

0 072 

- 

0-811 


0-003 

- 

- 

0-005 

0 003 

0 001 

- 

0-001 

0-987 


* Key: 1 financial, real estate, and insurance; 2 professional services; 3 commercial residential; 

4 general retail; 5 personal and commercial services; 6 eating and/or drinking; 7 parking; 

8 entertainment; 9 other. 





= known State Vector at T. 






Estimated second-order probability transition matrices of commercial activities for downtown Denver. First row, 1947-1951 -1963; second 
row, 1951-1963-1967; third row, 1963-1967-1971. Key as in table I. 

Commercial residential. General retail. 


Transition lo 

Transition to (second period) 
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(first period) 

1 

2 

3 

4 

5 

6 

7 

8 

9 

1 

2 

3 

4 

5 

6 

7 

8 

9 
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_ 

_ 

_ 

_ 


_ 

_ 

0*438 

_ 

_ 

0*281 

0*129 


0 031 

_ 

0*125 

real estate. 
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- 

- 

- 

- 

- 

- 

- 

0*750 

0*028 

- 

0*083 

- 


0*028 

- 

0*111 

and insurance 
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- 

- 

- 

- 

- 

- 

- 

0*889 

- 

- 
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0*055 



- 

- 

Professional 

. 



_ 

_ 

_ 

- 

_ 

_ 
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0*400 
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_ 

- 

services 
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- 

- 

- 

- 

- 

- 

- 

0*400 
1 00 

0*200 
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- 



- 

0*400 

Commercial 

0*051 


0*795 

_ 

m • 

. 

0*103 


0*051 


_ 

_ 





_ 


residential 

- 


1*00 

- 

- 


- 

- 

- 

- 

- 

1*00 


- 



- 

- 


- 


0*882 
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- 

- 
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- 

- 

- 

- 

- 


- 



- 

- 

General retail 

_ 


. 

_ 

_ 


_ 

_ 

1*00 

0 064 

0*010 

0*002 
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0*043 

0*030 
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0*002 
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- 

- 

- 

- 

- 

- 

- 
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0*003 
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- 

1*00 

- 

- 

- 

- 

- 
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_ 

• 

_ 

- 

_ 

_ 


_ 

_ 

_ 

0*080 

- 

0*080 
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_ 

0*240 

_ 

0*160 
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- 

- 

- 

- 

- 

- 

- 

- 

- 
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- 
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- 
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- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

0*200 
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- 
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- 


_ 

- 

_ 
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- 


_ 

_ 

_ 
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_ 
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. 

_ 

- 
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- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 
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- 
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- 
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- 

- 

- 

- 

- 

1*00 

- 

- 

- 

- 

- 

- 
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- 
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- 

- 
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- 

- 

- 

- 

- 

_ 

- 

- 

- 

_ 

- 

- 
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- 

. 
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- 

- 

- 

- 

- 

- 

1 *00 

- 

- 

- 

- 

- 
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- 
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- 

- 
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- 

- 

- 

- 

- 
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- 

- 
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- 
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- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

_ 

- 

- 

- 
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- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 
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- 
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0*238 0*020 0*009 0195 - 

0*167 - 0*009 0*238 0*026 
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0 432 
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l *00 


0*073 0*013 
0*101 0*004 
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DISCRIMINANT ANALYSIS 


A METHOD OF PRODUCING HYBRID VARIABLES SO AS TO PRODUCE THE BEST 
POSSIBLE SEPARATION, OR DISCRIMINATION, BETWEEN VARIOUS GROUPS. 
SUCH DISCRIMINANT FUNCTIONS ARE OFTEN NECESSARY BECAUSE OF CO- 
LLINEARITY AMONG INDEPENDENT VARIABLES ORIGINALLY USED IN THE 
ANALYSIS. 

THE METHOD INVOLVES TWO SETS OF EQUATIONS: (i) A SET RELATING 

THE GROUP MEMBERSHIP TO THE DISCRIMINANT FUNCTIONS; AND (ii ) A SET 
RELATING THE ORIGINAL VARIABLES TO THE DISCRIMINANT FUNCTIONS. 


1. TESTING (AND GENERATING) HYPOTHESES, e.g., INVESTIGATE WHETHER 
A CERTAIN SET OF RELATED VARIABLES (MEASURED ON INTERVAL OR 
RATIO SCALES) SUCCESSFULLY DISCRIMINATE BETWEEN GROUPS OF 
OBSERVATIONS ON A NOMINAL SCALE. 

2. EVALUATING A CLASSIFICATION, e.g., INDICATE THE NUMBER OF 
MISCLASSIFICATIONS IN THE DEPENDENT VARIABLE. 

3. ESTIMATING VALUES FOR OTHER OBSERVATIONS, i.e., ALLOCATING 


NEW OBSERVATIONS TO AN EXISTING CLASSIFICATION. 



Derivation of a discriminant function, showing: (A) the location of 12 
observations, in twogroups, on two orthogonal independent variables; (B) the location 
of members of the two groups, and the group means, on the separate independent 
variables; (C) location of a discriminant function which achieves maximal separation 
of the two groups; and (D) the location of members of the two groups, and the group 
means, on the discriminant function. 

L 




A discriminant function separating two groups of observations according to 
their values on four, related variables. 



A discriminant analysis showing the use of two discriminant functions to 
separate three groups of observations in a three-variable space. 
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